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Abstract. This paper deals with the solution of unified fractional reaction- 
diffusion systems. The results are obtained in compact and elegant forms in 
terms of Mittag-Leffler functions and generalized Mittag-LefHer functions, which 
are suitable for numerical computation. On account of the most general char- 
acter of the derived results, numerous results on fractional reaction, fractional 
diffusion, and fractional reaction-diffusion problems scattered in the literature, 
including the recently derived results by the authors for reaction-diffusion mod- 
els, follow as special cases. 

1 Introduction 

Reaction-diffusion models have found numerous applications in pattern forma- 
tion in biology, chemistry, and physics.. These systems show that diffusion can 
produce the spontaneous formation of spatio-temporal patterns. A piecewise 
linear approach in connection with the diffusive processes has been developed 
by Strier et al. (1995). A similar approach is recently used by Manne et al. 
(2000) to investigate effects on the propagation of nonlinear wave fronts. 

The simplest reaction-diffusion model can be described by an equation 



where D is the diffusion constant and F(N) is a nonlinear function representing 
reaction kinetics. It is interesting to observe that for F{N) — '-fN{l — N), (1) 




+ XF{N), N = N{x,t) 



(1) 
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reduces to the Fisher-Kolmogorov equation and if, we set F{N) = 7iV(l — A''^), 
it reduces to the real Ginsburg-Landau equation. 

A generalization of (1) has been considered by Manne et al. (2000) in the 
form 

d^N (IN ^d^N , 

where ^ indicates the strength of the nonlinearity of the system. 

In this article, we present a straightforward method for the systematic deriva- 
tion of the solution of nonlinear reaction-diffusion equations connected with non- 
linear waves, which is a generalization of the eq. (2). The results are derived 
in a closed form by the application of Laplace and Fourier transforms, which 
are suitable for numerical computation. The present study is a continuation of 
our investigations reported earlier in the articles of Saxena et al. (2002, 2004, 
2004a, 2004b, 2006a, 2006b, 2006c). 



2 Auxiliary Results 

A generalization of the Mittag-LefHer (1903, 1905) function 

^^(^) = E rrna + lV ("eC,i?e(a)>0) (3) 
was introduced by Wiman (1905) in the general form 

oo „ 

EaA^) = J2 w e C,Re{a) > 0,Re{P) > 0). (4) 

The main results of these functions are available in the handbook of Erdclyi, 
et al. (1955, Section 18.1) and the monographs by Dzherbashyan (1966, 1993). 
Prabhakar (1971) introduced a generalization of (4) in the form 

^^A'^ = T. j[na + (i){n)V 7 € C; iie(a), i?e(/3), i?e(7) > 0), (5) 

where (7)^ is the Pochhammer symbol, defined by 

(7)0 = 1, (7)r = 7(7 + 1)(7 + 2) . . . (7 + r - 1), r = 1, 2, . . . , 7 7^ 0. 

It is an entire function with p — [i?e(Q;)]~^ (Prabhakar, 1971). Solution of 
generalized Volterra-type differintegral equations associated with this function 
as a kernel is derived by Kilbas et al. (2002). A general theory of generalized 
fractional calculus based on this function has been developed by Kilbas et al. 
(2004), generalizing the results for Riemann-Liouville fractional integrals and 
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derivatives, which form the backbone of fractional differential equations. 
For 7 = 1, this function coincides with (4), while for /? = 7 = 1 with (3): 

El^ = E^,l3{z), El^{z)=E^{z). (6) 

We also have 

$(a,/3;^) =1 F^{a;j3;z)=V{(})Elp{z), (7) 

where $(a, /3; z) is Kummer's confluent hypergeometric function defined in 
Erdelyi et al. (1953, p. 248, eq. (1)). Prabhakar (1971, p. 8, eq. (2.5)) has 

shown that 

L {p-^e-''El^{ujt^)-s] = - LJS-^)-\ (8) 



where i?e(/3) > 0, i?e(7) > 0, Re{s) > and s > \uj\ "^w^) . 

In terms of the H-function, (5) can be expressed in the form 



KA-^) = fM^ii Hi7uU.)\ . (9) 

where Re{a) > 0, i?e(/3) > 0, Re{j) > 0. 

We will need the definitions of the Laplace and Fourier transforms of a 
function N[x,t) and their inverses, which are given below; 
The Laplace transform of a function N{x, t) with respect to t is defined by 



(10) 



roc 

L{N{x,t)} = N~{x,s) = e-'*N{x,t)dt, {t>Q),{xeR), 

Jo 

where Re{s) > 0, and its inverse transform with respect to s is given by 

L-^{N~{x,s)} = L-^{N^{x,s);t}=— e'*N^{x,s)ds, (11) 

7 being a fixed real number. 

The Fourier transform of a function N{x, t) with respect to x is defined by 

/oo 
e'^''N{x,t)dx. (12) 
-00 

The inverse Fourier transform with respect to k is given by the formula 

1 

F-^{F*{x,t)} = — e-''''=F*{k,t)dk. (13) 

The space of functions for which the transforms defined by (10) and (12) exist 
are denoted by LF = L{R+) x F{R). In view of the results (Mathai and Saxena, 
1978, p. 49), also see (Prudnikov et al., 1989, p. 355, eq. 2.25.3.2), the cosine 
transform of the H-function is given by 
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pOO 

Jo 



at" 



dt 



^ Trn+l,m 
^-"g+l,p+2 



where 



fc^ (l-fa„B,),(i+£,f) 
a (/,,/*),(l-ap,Ap),(i+£,f) 

I mas; 



(14) 



< 



> 1; fc, M > 0; i?e [p + /i™<7<„ (^) 

\arg a\ < ^n0; = " + ^7=1 " EI=„^+/s. > 0- 

The Riemann-Liouville fractional integral of order v is defined by (Miller 
and Ross, 1993, p. 45) 



oDi''N{x,t) 



n^) Jo 



{t - u)''-^N{x,u)du, 



(15) 



whore Re(h') > 0. 

Following Samko et al. (1990, p. 37), we define the fractional derivative of 
order a > in the form 



oD?N{x,t) = 



1 /■* 
i-a)dt^ Jo 



N{x,u)du 

-t-^^-"; - Y{n-a)dt^Jo (i -«)«-"+! ^ 



(t > 0), in = [a] + 1), (16) 



where [a] means the integral part of the number a. 
Prom Erdelyi et al. (1954b, p. 182), we have 



(17) 



where Re{s) > and Re(a.) > 0. 

The Laplace transform of the fractional derivative, defined by (16), is given 
by Oldham and Spanier, (1974, p. 134, eq. (8.1.3)) 

n 

L{oD^N{x,t)} = s"'N~{x,s) - oDt'''N{x,t)\t=o, n - 1< a < n. 

r=l 

(18) 

In certain boundary-value problems arising in the theory of viscoelasticity and 

in the hereditary solid mechanics, the following fractional derivative of order 
a > is introduced by Caputo (1969) in the form 



oD?f{x,t) 



1 



f^"'\x,T)di 



T{m - a) Jo (t - t)"+i-™ 

TO — 1 < a < m, Re{a) > 0, to e N 



if ( 



(19) 
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Caputo (1969) has given the Laplace transform of the derivative as 



m — 1 



L{oD«/(a;,i)} = s«r(a;,s)- ^s"-'-V^''^(a;,0), (m - 1< a < m), (20) 

i-=0 

where /~(.x, s) is the Laplace transform of f{x,t). 

The above formula is useful in deriving the solution of differintegral equations 
of fractional order governing certain physical problems of reaction and diffusion. 

We also need the Weyl fractional operator, defined by 

where n = [/i] + 1 is an integral part of /i > 0. 

The Fourier transform of the Caputo derivative is given by (Metzler and 
Klafter, 2000, p. 59, A.ll) 

F{_^D!if{x,t)} = {ikrnk,t). (22) 

Following the convention initiated by Comptc (1996), we suppress the imaginary 
unit in Fourier space by adopting a slightly modified form of the result (21) in 
our investigations (Metzler and Klafter, 2000, p. 59, A. 12): 

F{_^D^;:f{x,t)} = -\k\^r{k,t). (23) 

3 Inverse Laplace Transform of Certain 
Algebraic Functions 

This section deals with the evaluation of the inverse Laplace transforms of cer- 
tain algebraic functions which are directly applicable in the analysis that follows. 
It will be shown here that 

IT—. — ;tUy(-l)'- 

^ r=0 

X:([)a'6'-'i(«-^)'-+(--^)'+«--^;+i_^),+(^_,),+„_^+,(-ct«), (24) 



1=0 

where 



Re{a) > 0,Re{P) > 0,Re{j) > 0,Re{s) > 0,Re{p) > 0, ^^^^ < 1, E'g J.) 



is the generalized Mittag-LefHer function defined by eq. (5) and provided that 
the series in eq. (24) is convergent. 
Proof. Assume that a > 7 > /3. We have 
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s" + asl^ + 6sT + c 



oo 

S 

r=0 

°° „p+(7-Q)r+(/3-7)i-a-l 



„p-^-iy> (-ir(a + 6a^-/^r 
^ (s«-/3 + cs-/5)'-+i 



= E(-irE(FV^''"' + • (25) 

r=0 !=0 \ - J 

On term by term inverting (25) with the help of the formula (8), it readily gives 
the desired result. The term by term inversion is justified by virtue of a theorem 
due to Doetsch (1956, §22). 

(i) From (25), it readily follows that 

L-i <^ — --t \ = e-f'Y{-aYt^'^~^>E^+\, ^. (-6t"), (26) 

I _|_ „ „/3 _|_ J, ' I /.^^ a,OL+(a-P)r-p+l\ /' V ; 



s°+6 



< 1 and provided 



where Re{a) > 0, Re{f3) > 0, Re{p) > 0, Re{s) > 0, 
that the scries in (26) is convergent. 

Some special cases of the result (26) are worth mentioning. 

For p = a, (26) reduces to the following result given by the authors (Saxena 
et al., 2006b): 

f „a-l ^ °° 

s^+as^+b ''\ = T.(--rt'"-'''^Zii-,y+.i-bn, (27) 

' ' r=0 

where Re{a) > 0, Re{(3) > 0, Re{s) > 0, \ ^^\ < 1 and provided that the series 
in (27) is convergent. 

Further, if we set p = P, then we arrive at another result given by the 
authors: 

L-i \ —-5 — -]t )■ = y (-a)¥"-'^)''i;''y (-6^), (28) 

I + /js/3 + ' I Z-J^ ' a,(a-/3)(i-+l)+lV n \ I 

^ ' r=0 

< l,a > /?, and provided that 



where Re{p.) > 0,i?e(/3) > 0,Re{s) > 0, 
the series in (28) is convergent. 



(ii) If we set /> = 1 in (24), we obtain the following result: 



\ s"' + asl^ + bs^ + c' J 

^ r=0 
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'y^^^-j^;^r-i^(a-7)r+(7-/3)i+0!-l 



X E 





r+1 



a,a+(Q— /3)r+(7— /3)/+a 



(-Ci"), 



where i?e(a) > 0, Re{(3) > 0, i?e(7) > 0, Re{s) > 0, 
that the series in (29) is convergent. 



(29) 

< 1, and provided 



(iii) Next, if we set p = a in (24), we then arrive at the following result 

oC« — 1 



1=0 



where i?e(a) > 0, Re{f3) > 0, i?e(7) > 0, Re{s) > 

that the series in (30) is convergent. 

(iv) Similarly, p = (3 in (24) will yield the formula 



< 1, and provided 



-\t\ 



3°" + asf^ + bs^ + c 

r 

(31) 

< 1, and provided 



r=0 



1=0 



where Re{a) > 0,i?e(/3) > 0,i?e(7) > 0,Re{s) > 0, 
that the series in (31) is convergent. 

(v) When p = 7,(24) readily gives the formula 



L-1 



^7-1 



+ as^* + bs'f + c 

= E(-l)^E(F)«'''''"'*^""'^''^^'"'^'^""'^«^a-7).+(7-/3)^+a-7+l(-^*")' 

(32) 

< 1, and provided 



1=0 



where Re{a) > 0, i?e(/3) > 0, i?e(7) > 0, i?e(s) > 0, 

that the series in (32) is convergent. 

(vi) When p = a + (3 then (24) reduces to 



ga+/3-l 



s« + as^ + 6pT + c 



^4 
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r=0 1=0 



(33) 
< 1, and 



where Re{a) > 0, Re{P) > 0, Re{j) > 0, Re{a + /?) > 
provided that the series in (33) is convergent. 



4 A General Case 

In this Section, it will be shown that if Re{aj) > (j = 1, . . . , n), Re{s) > 0, 
then 



1 



m! 



n 



ri >0...rn>0 



E"=i{(°'-+i/°i)''°''+'"°'} 

(s<*i-<»2 + (ao/ai)s-02) 



where A = J2"^-^{a2-aj+i)rj+ai{l+m)-a2m, 

provided that the scries in (34) is convergent. 

Proof: Let us assume that a2 > cxj+i, (j = 2, 3, . . . , n). We have 

sj^ 

ao + ais«i + a2s"2 + ass"^ + . . . + a„s"" + a„+is"+i 



(34) 
< 1, 



[s«i-«2 + (ao/ai)s-"2 

oo 

-1^^-02-1 (— 1)™ 

m=0 



1 I ELJ(°-+i/°i)^°'-+^"°^} 



(s°l-°2+(oo/ai)s-°2) 

(sai-a2 + (ao/ai)s-"2) 



m+1 ■ 



(35) 



On applying the multinomial theorem (Abramowitz and Stegun, 1968, p. 823, 
para 24.1.2): 

{Xi + X2 + ■ ■ ■ + X,,,)" = E("!"l'"2, • ■ . ,nm)x"'X2^ . . . X^"" , 

summed over rii + n2 + ■ ■ ■ + rim = ri; the above line transforms into the form 



:^E(-ir E 



ai — ' ■ — ' ri ! . . . r~ M 

m=0 ri + ...+r„ = m 

ri>0,...,rn>0 



TT" 



ai 
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y^,"_j (Q3 + l-a2)rj+p+a2rra-(l+"t)o!i-l 



(1 + (ao/ai)s-«i) 



m+l 



(36) 



Interpreting the above expression with the help of the formula (8), we obtain 
the desired result. This completes the Proof of (34). Some interesting special 
cases of the general result (34) arc given below. More results can be derived 
by specializing the parameters and the variable and applying certain known 
theorems of Laplace and inverse transforms, 
(i) If we set p = ai , then (34) reduces to 



L- 



+ a„s"" + a„+is«n+i 



X E 



m=0 



[Tn+l 

ai,A — ai + 1 



T->O...T-„>0 



ri!...r„! 



n 



ai 



i-A—ai 



ai 



(37) 



where 

i?e(aj) > (j = l,...,n),i?e(s) > 0, 



(s°l-°2+(ao/ai)s-»2) 



< 1,^ = 



~ + Q;i(l + to) — a;2TO, and provided that the series in (37) is 

convergent. 

(ii) On the other hand, if we take p = a2 in (34), it gives the formula 



ao + ais"i + a2S"2 + a„_2s"-^ + . . . + a„s"" + a„+is""+i 



TO! 



m=0 



y Jpm+1 

^ ^ai,A-a2+l 



T-i>0...r„>0 



n! . . .r„! 



TTTl 



ai 



iA-a2 



(38) 



where Re{aj) > (j = 1, . . . , n), i?,e(s) > 0, 



E"=i{(°'-+i/°i)"°''+^'°'} 

(5ai-o2+(ao/ai).s-<'2) 



< 1, 



A = X]J=i(q^2 — oij+{)rj + ai(l + to) — Q!2TO,, and provided that the series in 
(38) is convergent . 

(iii) Next, if we set p = 1 in (34), then we arrive at the following inversion 
formula given earlier by Podlubny (1999) in a slightly different form: 



1 



ao + ais"i + a2s"2 + a„_2s"3 + . . . + a„s"" + an+is^^+i ' 
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Oil, A 



ai 



(39) 



where 

Re{aj) > 0(j = 1, . . .,n),Re{s) > 0, 



E:=i{K+i/ai)«"'-+^-"n 



(s«i-«2 + (ao/ai)s-°2) 



< 1, 



A = X]J=i('^2 — + Q!i(l + TTi) — a!2?Ti, and provided that the series in 

(4.5) is convergent. 

(iv) When p = ai + a2, (34), then it reduces to the formula 

gai+a2-l 

■,t 

ao + ais«i + a2S°2 + a„_2S«3 + . . . + a„s«n + a„+is'^^+^ 



m=0 



ai,Ai+l 



ri+r2 + ---+r- 
r>0...r-„>0 



ri! . . . r„! 



n 



ai 



ai 



(40) 



where Re{aj) > 0.(j = 1, 2, . . . ,n) : Q2 > Q:j+i(i = 2, 3, . . . , n), ai 9^ 0, J?e(s) > 
0, 



E"=i{K+i/«i)s"-^-"^} 



(s«i-«2 + (ao/ai)s-«2) 



< 1, Ai = E(^2 - Q!j+i''j + aim - a2(TO + 1), 



and provided that the series in (40) is convergent. In what foUows, ,.^(.) will 
be employed to denote the generalized Mittag-Leffler function, defined lay (5). 



5 Solution of Fractional Reaction-Diffusion 
Equation 

In this Section, it is proposed to derive the solution of the fractional diffusion 
system connected with nonlinear waves governed by the eq. (41). This system is 
a generalized form of the reaction-diffusion equation recently studied by Manne 
et al. (2000). The result is given in the form of the following 
Theorem 1. Consider the fractional reaction-diffusion equation associated with 
Caputo derivatives in the form 

oDl'N{x, t) + ao DfN{x, t) + boD]N{x, t) 
= -ooD2N{x,t)+fN{x,t)+(P{x,t) (41) 
(0<a<l,0</3<l,0<7<l,r?>0) 
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with initial conditions 

N{x, 0) = f{x), xe^, lim^^±ooN{x, t)=0,t>0, (42) 

where u"^ is a diffusion constant, ^ is a constant which describes the nonhncarity 
in the system, and <f{x,t) is nonhnear function which belongs to the area of 
reaction-diffusion, then there holds the following formula for the solution of 
(41): 

N{x,t) = 7rl2^-^y / 6"^'=" V([)a'6'-'t("-'')''+(^-'^)V*(A:) 

{-^Q^a-7)r+(7-/3)i+l(~'^*") 

' Q:,(a — 7)r+(7 — p)/ + lV / Q>^(^a—^)r-\-['y—f3)l-\-l\ /J 

(43) 

where c = — and provided that the series and integrals in (43) are 

convergent. 

Proof. Applying the Laplace transform with respect to the time variable t and 
using the boundary conditions, we find that 

s"iV~(.T, .s) - s"-\f (x) + as'^N'-ix, s) - as^-^f{x) + 
bs^N~{x,s) -bs'^-^fix) 
= ly'^ -^U;iN-'{x,s)+^^N-{x,s)+^~{x,s) (44) 

If we apply the Fourier transform with respect to the space variable x, it yields 

s"N~*{k, s) - 3°'-^f*{k) + asf^N'-*{k, s), -as^-'^f*{k) + 
bs'^N~*{k,s) - 6s^-V*(fc) 
= -iy^\k\'^N*{k,s)+^^N~*{k,s) + 'fi^*{k,s). (45) 

Solving for N'"*{k, s), it gives 

N-*(ks)= i^"-'+<^^'-'+bs--')r{k)+^-*{k,s) 

^ ' ' 8°^ + asf^ + bs'i + c ^ ' 

where c = h''^\k\^ — 

To invert the eq. (46), it is convenient to first invert the Laplace transform 
and then the Fourier transform. Inverting the Laplace transform with the help 
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of the results (29), (30), (31) and (32), it yields 

OO r 

N*{k,t) = ^(-l)'-^([)a'6'-'t(«-T)''+(T-'5)'/*(A;) 



r=0 



1=0 



fpr+1 ( _r.+°' \ 

Q,(a-7)r+(7-/3);+lV " > 

,a-l3 fpr+l 



+ at" -B^ („_^)^_,_(^_^);_,_„_^_,_i( ct") + 
bt" '''-E^^^_^)r_,_(^_^);_,_„_^_,_i(-Ci") 



r=0 ;=o 

/* ^*{k, t - o^^"-^^''+(^-''^'+"-^i^;:^^_,).+(,_^),+„(-cr)d6 

J 



(47) 



Finally, the inverse Fourier transform of (47) gives the desired solution in the 
form 

, oo .oo r 

N{x,t) = — Y.i-'^y / e-''^^ y^{l)a'b'-¥'^-^'>'~+^^-f^'>'f*{k) 



r=0 

^a,{a-j)r+{-y-0)k+l\ " ) 



dk 



a,(a— 7)'"+(7— /3)i + l ^ ' 



r=0 ^'^ ;=o "^0 

/oo 
e-^'=>*(fc,t-e)£;;+;^_^),^(^_,),^„(-cr)dA, 
-OO 

This completes the proof of the theorem. 



(48) 



6 Interesting Special Cases and Fundamental 
Solutions 

Some special cases of the theorem 1 are given below, (i) If we set a = (3 = •y ■ 
rj — 1/2, the theorem reduces to 

Corollary 1.1. Consider the fractional reaction-diffusion equation 



oD]'^N{x, t) + ao d]'^N{x, t) + boD]'^N{x, t) 
= ,ylDl/^N{x, t) + eN{x, t) + <p(x, t), 



(49) 
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with initial conditions 

N{x, 0) = f{x),x e lim:c^±ooN{x, t)=0,t>0, (50) 

where is a diffusion constant, ^ is a constant which describes the nonhncarity 
in the system, and is nonhnear function which belongs to the area of reaction- 
diffusion, then there holds the following formula for the solution of (49): 

N{x,t) = (l + a + ^) J2i-iy r e-^>'-j2iiWb'--'r{k)El%\{-c*ndk 

^'^ r=a 1=0 

r-=o ^'^ ;=o "^0 
X El+l^^^{-c*^^)dkd^ (51) 

where c* = v'^\k\^''^ - . 

(ii) When f{x) = S{x), where S{x) is the Dirac-delta function, the theorem 

yields the following result 

Corollary 1.2. Consider the fractional reaction-diffusion system governed by 
the equation 

oD^N{x, t) + aoD^N{x, t) + boD]N{x, t) 
= ul^D2N{x,t)+eN{x,t)+cpix,t) (52) 
(0 < a < 1,0 < /3 < 1,0 < 7 < 1,7? > 0), 

subject to the initial conditions 

N{x, 0) = 5{x),x G 3?, lim^^±ooN{x, t) = 0,t > 0, (53) 

where 6{x) is the Dirac-delta function. Here ^ is a constant that describes the 
nonlinearity in the system, and (p{x,t) is nonlinear function which belongs to 
the area of reaction-diffusion. Then there exists the following formula for the 
fundamental solution of (52) subject to the initial conditions (53): 

r=0 1=0 

a,{a—^)r+{-y—/3)k+ll^ ' 
^a^i-7)r-h(7-/3);-Hl(~'^*")} 



r=0 1=0 "^0 
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f 

o — ( 



^;^L,).+(,_^),+„(-c^")rffcrfC (54) 

where h = v'^\k\^ ~ and provided that the series and integrals in (54) are 
convergent. 

As 6 ^ 0, we obtain the following result recently obtained by the authors 
(Saxena et al., 2006b): 

Corollary 1.3. Consider the fractional reaction-diffusion equation 

oD^N{x,t)+aoD^tN{x,t) 
= ul^D2N{x,t)+eN{x,t) + ^{x,t) (55) 
(0 < a < 1,0 < /? < 1,7 > 0) 

with initial conditions 

N{x, 0) = /(x), e 3?, lim^^±^N{x, t)=0,t>0, (56) 

where is a diffusion constant, ^ is a constant which describes the nonlinearity 
in the system, and (j){x,t) is nonlinear function which belongs to the area of 
reaction-diffusion, then there holds the following formula for the solution of 
(55). 

N{x,t) 

„ n J—oo 

dk 



r=0 



+1 

a+(a— /?)r 



{-bC)dkd^, (57) 

where a > /3 and b = v^lkp —S,"^. Now if we set f{x) = S{x), = 2, a is replaced 
by 2a and /3 by a, 7 = 0, and use the result (Saxena et al., 2006b, eqs. (24) and 

(28)), the following result is obtained. 

Corollary 1.4 Consider the following reaction-diffusion system 
d'^°'N(x,t) d°'N(x,t) ^d'^N{x,t) x / , r. ^ ^ ,.c^ 

with the initial conditions 

N{x,Qi) = 5{x),x£^,Nt{x,Qi) = Q, lim N{x,t) = Q,t > Q, (59) 

cc— '■ztoo 



14 



4>{x, t) is a nonlinear function belonging to the area of reaction-diffusion. Then 
for the fundamental solution of (58) subject to the initial conditions (59), there 
holds the formula 



1 



+ 00 



N(x,t) = i[ / exp(—ikx) 

^ 27r 7(^^^46)7-00 

{(A + ojEaiXt") -{fJ. + a)Ea{fJ.t°')}dk 

1 pt r + co 

+ expi-ikx)v*{k,t-0 

X [S„,„(Ar) - E^,^{fiC)]dkd^, (60) 
where A and n are the real and distinct roots of the quadratic equation 

y^ + ay + b = 0, (61) 

given by 

A = li-a + V(a2 - 46))and = ^{-a - ^(a^ - 46)), (62) 

where 6^ = i^^fc^ — and provided that the integral appearing in (60) is con- 
vergent. 

Next, if wc set (p{x, t) = 6 = 0, 77 = 2, replace a by 2a, and /3 by a in (58), we 
then obtain the following interesting result, which includes many known results 
on fractional telegraph equations including the one recently given by Orsingher 
and Beghin (2004). 

Corollary 1.5. Consider the following reaction-diffusion system with the initial 
conditions 

d^"N(x,t) d"N(x,t) ^d'^N(x,t) ,2.r/ x ^ 

with the initial conditions 

N{x,0) = 6{x),x G^,Nt{x,0) =0 lim N{x,t) = 0,t > 0, (64) 

Then for the fundamental solution of (63) subject to the initial conditions (64), 
there holds the formula 

1 /■+°° 
N{x,t) = - — ====[/ exp{-ikx) 
2-K^{a^ - 46) i-oo 

{(A + a)E^{X') - (/X a)Eo,{iit°')\ dk, (65) 

where A and are defined in (62), 6 = v^k"^ —^^ and Eoi{x) is the Mittag-LefHer 
function defined by (3). 
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If we set = 0, then the Corollary 1.5 reduces to the result, which states 
that the reaction-diffusion system (66) 

d'^"N{x,t) d°'N(x,t) ^d^N{x,t) „ 

^ o + Q ^ ^ =u' — ^V^, 0<a<l (66) 

with the initial conditions 

N(x,t) = 5{x),x & U,Nt(x,0) = Q, lim:o^±^N{x,t) = 0,t>0, (67) 
has the fundamental solution, given by 



1 



+ 00 



{{X + a)Ea{Xt°') - {fi + a)Ea{ne)}dk, (68) 

where A and fi are defined in (62), b = v'^k'^ and Ea{x) is the Mittag-Leffler 
function defined by (3). 

The result (68) can be rewritten in the form 

1 /■+°° 
N{x,t) = — / exp{—ikx) 

where A and fj, are defined in (62) and is the Mittag-Leffler function, defined by 
(3). 

The eq. (69) represents the solution of time-fractional telegraph eq. (66) 

subject to the initial conditions in (67), recently solved by Orsinghcr and Bcghin 
(2004). It may be remarked here that the solution as given by Orsingher and 
Beghin (2004) is in terms of the Fourier transform of the solution in the form 
given below. It is interesting to observe that the Fourier transform of the solu- 
tion of the eqs. (66) and (67) can be expressed in the form 

N*{x,t) = U{i + ^=^==)E^{\n + (1 - ^^^=)E^{^,n\ , 

(70) 

where A and /x are defined in (62) and Ea{x) is the Mittag-Leffler function 
defined by (3). 

Following a similar procedure, the following general theorem can be estab- 
lished which makes use of the general result (24). 
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7 General Theorem 2 

Consider the unified fractional reaction-difi^usion equation associated with Ca- 
puto derivatives 

ai oD^'N{x, t) + 02 qD'^^N{x, t) + . . . + a„+i oD'i"^^N{x, t) 
= u"^ -ooD^N{x,t)+^^N{x,t) + (pix,t) (71) 
< < 1, (i = 1, . . . , n + 1), > 

with initial conditions 

N{x, 0) = /(x), a; e 5R, ;im^^±ooA^(a;, = 0, t>0 (72) 

where is a difl:usion constant, ^ is a constant which describes the nonlinearity 
in the system, and (f){x,t) is nonlinear function which belongs to the area of 
reaction- diffusion, then there holds the following formula for the solution of 
(71). 

N{x,t) = -—^(-1)-/ r{k)exp{-ikx) Yl 

^'^"'1 m=0 ri+.. .+r„=m 



n! 

n+l 



aj=l 

+ ^ E(-ir/ / ea;p(-^Me^-V*(fc,<-e) 



ri +...+r,, — m 



ai 



X K;^li~t"')dkd^, (73) 

where ao = i''^\k\^ — ^'^,Re{aj) > 0,)j = 1, . . . ,n);Re{aj+i — a2) > 

for j = 2, . . . ,n,limx^±coN{x,t) = 0; A = J2'j=i{(^2-aj+i)rj+ai{l+m)-a2m 
and provided that the series and integrals in (73) are convergent. 

8 Consequences of General Theorem 2 

When aj = iJL = l/2{j = 1, 2, . . . , n = 1), the theorem 2 yields 
Corollary 2.1. Consider the unified fractional reaction-diffusion equation as- 
sociated with Caputo derivatives 

ai oDI^^N{x, t) + a2 oDI'^N{x, t) + ... + a„+i oD]'^N{x, t) 
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= -^Dl/^N{x, t) + e'iV(x, t) + ^{x, t) (74) 
0<aj<l,0- = l,...,n+l) 

with initial conditions 

N{x, 0) = f{x), {x e 3?), lim^^±ooN{x, t)=0,t>0 (75) 

where is a difTusion constant, ^ is a constant which describes the nonhncarity 
in the system, and (j){x,t) is nonUnear function which belongs to the area of 
reaction-diffusion, then there holds the following formula for the solution of 
(74) and (75). 



N{x,t) 



1 oo ^oo 

^ m=0 -'-oo ri+...+r„=m ^ 



n}Li ( ^ ) ' } (ai + . . . + an+i)E^}+^^{-^)dk 



m=0 

Irrn ( 

ri! 



where b = v'^\k\^/'^ — ^^,ai ^ 0, and provided that the series and integrals in 
eq. (76) are convergent. 

Corollary 2.2. Consider the fractional reaction-diffusion equation 

ai oD^'Nix, t) + 02 qD-^'N{x, t) + . . . + a„+i qD'^"+^N{x, t) 
= 1/2 _oo^^iV(x,f)+^2 ^(x,t)+v?(x,i) (77) 
< ttj < 1, (j = 1, . . . , n + 1), > 

with initial conditions 

N{x, 0) = 5{x), (a; e K), lim^^±ooN{x, t) =0,t>0, (78) 

where is a diffusion constant, ^ is a constant which describes the nonlincarity 
in the system, and '^{x,t) is nonlinear function which belongs to the area of 
reaction-diffusion, then there holds the following formula for the fundamental 
solution of eqs. (77) and (78). 



1 °° /-oo 



m=0 ri+rn=m 



n ! . . . n 
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til / J -, f^i 



1 ^ roc 



ri + ...+r„=m v. \ 



X (79) 

where A = X]J=i(Q^i+i - 0^2)?"^ + cci + (ai - 02)771, Re{aj) > 0, (j = 1, . . . , n); 
i?e(Q:j+i — 02) > 

for j = 2, . . . n, ao = i'^|A:|'^— lima;^±oo-^(a;, t) = 0, ai ^ 0, and provided 
that the integrals and series in (79) are convergent. 
On the other hand, if wc set ^ = 0, wc arrive at 
Corollary 2.3. Consider the fractional reaction-diffusion equation 



ai oD^'Nix, t) + a2 o^t"' ^(a^, *) + ••• + a^+i A""^' ^(^: 
= -ooD^N{x, t) + ^{x, t) (80) 
0<Q!<l,0</3<l,0<7<l,/x>0 

with initial conditions 

N{x, 0) = 5{x), {x e SR), Zim^^±ooiV(^, t) = 0, i > 0, (81) 

where i^^ is a diffusion constant, and 4'{x, t) is nonlinear function which belongs 
to the area of reaction-diffusion, then there holds the following formula for the 
fundamental solution of eqs. (80) and (81). 

m=0 n + .-.+rn^m k \ / j 

Va^ / Gi(x-r,<)/(T)dT 

m=0 " ri4-...-|-r„=m ^ \ / ^ 

X r G2{x-T,t-0^{T,0dTd^, (82) 
JO 

where, A = - Oi2)rj + ai + (ai - oi2)m, Re{aj) > 0, (j = 1, . . . , n); 

iie(Q!j+i — 02) > 
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for j = 2, . . . ,n,ao = i'^\k\'^ — limx^±<xN{x,t) = 0, and ai ^ 0; 



Gi{x,t) 



'X— >±CXI-i 

im+1 



ai 



dk 



— - — / cos(kx)Hl'2 
1/m 



( — m,l/ fi 

(0,l//.),(c«.-^,^) 



■A—au 



t 

/i(m!)|a;| 



■2,1 



3,3 



a; o 



+ (1,1/2) 
(1,1), (l+m,l/,*), (1,1/2) 



(83) 



and 



Gi{x,t) = exp{-ikx)E:^+l 



yn(TO)! 



cos{kx)H\'2 



ai 

j,2/Miai/M|;!j| 



(-m,l//t) 

(0,l/;t*),(l-A,c<i/,i) 



/i(m)!|a; 



2,1 
3,3 



l/Ml 



1/2 //i 



(l,l/M),(A,ai/M), (1,1/2) 
(1,1), (l+m,l/,i), (1,1/2) 



(84) 



provided that the series and integral in eq. (82) arc convergent. 



Remark. It is interesting to observe that the method employed for deriving 
the solution of the equations (41) and (42) as well as (71) and (72) in the space 
= LF = L{R^) X F{R) can also be apphed in the space LF' = L'{R^) x F' , 
where F' = F'(r) is the space of Fourier transform of generalized functions .As 
an illustration, we can choose F' = S' or F' = D'. The Fourier transforms 
in S' and D' arc introduced by Schwartz and Gclfand and Shilov, respectively. 
S' is the dual of the space S, which is the space of all infinitely difTcrcntiable 
functions which, together with their derivatives approach zero more rapidly than 
any power of l/\x\ as \x\ oo (Gclfand and Shilov. 1964, p. 16). Here D' is 
the dual of the space D which consists of all smooth functions with compact 
supports (Brychkov and Prudnikov, 1989, p. 3). For further details, the reader is 
referred to the monographs written by Gclfand and Shilov (1964) and Brychkov 
and Prudnikov et al. (1989), if we replace the Laplace and Fourier transforms 
in eqs. (10) and (12) by the corresponding Laplace and Fourier transform of 
the generalized functions. 
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